THE SANDPILE GROUP OF A TREE 



LIONEL LEVINE 

Abstract. A wired tree is a graph obtained from a tree by collapsing the 
leaves to a single vertex. We describe a pair of short exact sequences relat- 
ing the sandpile group of a wired tree to the sandpile groups of its principal 
subtrees. In the case of a regular tree these sequences split, enabling us to com- 
pute the full decomposition of the sandpile group as a product of cyclic groups. 
This resolves in the affirmative a conjecture of E. Toumpakari concerning the 
ranks of the Sylow p-subgroups. 



1. Introduction 

We begin with a simple combinatorial problem. Fix integers n > 2 and d > 3. 
By the d-regular tree of height n we will mean the finite rooted tree in which each 
non-leaf vertex has d — 1 children, and the path from each leaf to the root has 
n — 1 edges. We denote this tree by T n . The wired d-regular tree of height n is the 
multigraph T n obtained from T n by collapsing all the leaves of to a single vertex s, 
the sink, and adding an edge connecting the root r to the sink. We do not collapse 
edges; thus each neighbor of the sink except for r has d — 1 edges to the sink. The 
principal branches of T n are the subtrees rooted at the children of the root. 

Lemma 1.1. Let t n be the number of spanning trees ofT n . Then for n > 4, 

t n = t d n -_\(dt n . 1 -(d-l)t d n -_l). 
Proof. If the edge (r, s) from the root to the sink is included in the spanning tree, 
then each of the principal branches of T n may be assigned a spanning tree inde- 
pendently, so there are t d ~_\ such spanning trees. On the other hand, if (r, s) is 
not included in the spanning tree, there is a path r ~ X\ ~ . . . ~ x n _i = s in the 
spanning tree from the root to the sink. In this case, every principal branch except 
the one rooted at x\ may be assigned a spanning tree independently; within the 
branch rooted at xi, every subbranch except the one rooted at x-i may be assigned 
a spanning tree independently; and so on (see Figure 1). Since there are [d— l)™ -1 
possible paths X\ <~ . . . <~ x n _\, we conclude that 

n-1 

(1) tn=t d n Z\ + (d-lT- 1 l[t d -\ 

fe=l 

Substituting n — 1 for n we find that 

n-2 

(d-ir- 2 n*2- 2 = *n-l-t£2> 
fc=l 
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Figure 1 . The two cases in the proof of Lemma 



hence from ([I]) 



tn = tt-\{tn-l + (d- !)(*„-! - t d -_\)). □ 



From Lemma |1.1| one can readily show by induction that 

n-2 

(2) t n = (1 + a + . . . + a"" 1 ) Y[(l + a+... + a k ) a " 



'(a-1) 



fc=l 



where a — d — 1. A variant of this formula was found by E. Toumpakari |12j . who 
gives an algebraic proof. 

For any graph G there is an abelian group, the sandpile group, whose order is 
the number of spanning trees of G; its definition and properties are reviewed in 
section 2. A product formula such as Q immediately raises the question of an 
analogous factorization of the sandpile group. Our main result establishes such a 
factorization. Write Z^ for the group (Z/pZ) © ... © (Z/pZ) with q summands. 



Theorem 1.2. The sandpile group of the wired regular tree T n of 
and height n is given by 



d = a+ 1 



SP(f n ) ~ z a u 



\a-X) 



©Z" 



V-i) 



©Zi_ 



-a+...+a" 



'l+o+a 2 ^ ' ■ ' W ^l+a+...+a"- 

In [10] we give an application of this result to the rotor-router model on regular 
trees. 

Toumpakari [12] studied the sandpile group of the ball B n inside the infinite d- 
regular tree. Her setup differs slightly from ours in that there is no edge connecting 
the root to the sink. She found the rank, exponent, and order of the sandpile group 
SP(B n ) and conjectured a formula for the ranks of its Sylow p-subgroups. We use 
Theorem 1 1.2| to give a proof of her conjecture. 

We remark that Chen and Schedler [S] study the sandpile group of thick trees 
(i.e., graphs obtained from trees by replacing some edges with multiple edges) 
without collapsing the leaves to the sink. They obtain quite a different product 
formula in this setting. 
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The remainder of the paper is structured as follows. In section 2, we briefly 
review the definition and basic properties of the sandpile group and recurrent states. 
In section 3, we characterize the recurrent states on a wired tree explicitly in terms 



of what we call critical vertices. We prove a general result, Theorem 3.3 relating the 
sandpile group of an arbitrary wired tree T to the sandpile groups of its principal 
branches T 1 , . . . , T k . This result takes the form of a pair of short exact sequences 

(3) -> R(f) -> SP(f) -► H(f) -> 

k 

(4) R(f l ,...,f k ) -^(§SP(f l ) -> H(T) ->0. 

i=i 

The groups involved are defined in section 3. The addition of an edge from the root 
to the sink is crucial here: it plays the same role in the full tree that the edge from 
x to r plays in the branch rooted at x. 

In section 4, we show that the sequences Q and Q are split when T is a regular 
tree. This allows us to express the sandpile group SP(T n ) of the wired regular tree 
as the direct sum of a cyclic group and a quotient of the direct sum of d — 1 copies 



of SP(T n -i), which enables us to prove Theorem 1.2 by induction 



Finally, in section 5, we deduce Toumpakari's conjecture from our main results. 

2. The Sandpile Group 



Let G be a finite graph with vertices X\, . . . ,x n . We single out one vertex, x 
called the sink. The reduced Laplacian A of G is the ii-lxn-1 matrix 



n • 




where di is the degree of Xi, and dij is the number of edges connecting Xi and Xj 
(we allow multiple edges, but not loops). The sandpile group of G is defined as the 
quotient 

SP(G)=Z n ~ 1 /A Z"- 1 . 
This group was defined independently by Dhar [7] , motivated by the abelian sand- 
pile model of self-organized criticality in statistical physics [1J, and by Lorenzini 
PTj in connection with arithmetic geometry. In the combinatorics literature, other 
common names for this group are the critical group [3] and the Jacobian 0. 

The sandpile group can be understood combinatorially in terms of chip-firing 
[31 14]. A nonnegative vector u G Z n_1 may be thought of as a chip configuration 
on G with Ui chips at vertex Xi- A vertex x^ is unstable if Uj > dj. An unstable vertex 
may topple, sending one chip along each incident edge. Note that the operation of 
toppling the vertex Xi corresponds to adding the column vector Aj to u. We say 
that a chip configuration u is stable if every non-sink vertex has fewer chips than 
its degree, so that no vertex can topple. If u is not stable, one can show that 
by successively toppling unstable vertices, in finitely many steps we arrive at a 
stable configuration u° . Note that toppling one vertex may cause other vertices to 
become unstable, resulting in a cascade of topplings in which a given vertex may 
topple many times. The order in which topplings are performed does not affect the 
final configuration u°; this is the "abelian property" of abelian sandpiles [3 [8]. 

The operation (u, v) i— > (u + v)° gives the set of stable chip configurations the 
structure of a commutative monoid, of which the sandpile group is a subgroup. A 
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stable chip configuration u is called recurrent if there is a nonzero chip configuration 
v, such that (u+v)° = u. One can show that every equivalence class of Z™ _1 modulo 
A has a unique recurrent representative. Thus the sandpile group SP(G) may be 
thought of as the set of recurrent configurations under the operation (u, v) i— > 
(u + v)° of addition followed by stabilization. For proofs of these basic lemmas 
about recurrent configurations and the sandpile group, see, for example 0[9]. 
If v is a nonnegative configuration, its recurrent representative is given by 

v := (u + e)° 

where e is the identity element of SP(G) (the recurrent representative of 0); indeed, 
v is recurrent since e is recurrent, and v = v (mod A) since e G AZ" -1 . Note that 
if u is a recurrent configuration and v is a nonnegative configuration, then 

(5) (u + v)° = (u+(v + e)°)° = (ti + v + e)° = ((u + e)° + v)° = (u + v)° . 

We will need just one additional fact about recurrent configurations, a criterion 
known as the "burning algorithm" |7J that tests whether a configuration is recurrent. 
We include a proof for the sake of completeness. 

Burning algorithm. Let f3(x) be the number of edges in G from x to the sink. 
A stable chip configuration u on G is recurrent if and only if adding (3{x) chips at 
each vertex x causes every vertex to topple exactly once. 

Proof. Note that 

n-l 

(6) /3 = A n = -^Ai. 

i=l 

If every vertex topples exactly once in the stabilization of u + f3, then 

ra-l 

(u + 13)° = u + (3 + A » = u, 

i=l 

so u is recurrent. Conversely, suppose u is recurrent. Since (3 £ AZ™ _1 we have 
(3 = e, hence from ^ 

(u+/3)° = (u + e)° =u. 

By ([g]), since are linearly independent, every vertex topples exactly once 

in the stabilization of u + (3. □ 

3. General Trees 

Let T be a finite rooted tree, and let T be the graph obtained by collapsing all 
the leaves of T to a single vertex s, the sink, and adding an edge connecting the 
root to the sink. Denote by C(x) the set of children of a vertex x € T. We first 
characterize the recurrent configurations on T explicitly. The characterization uses 
the following recursive definition. 

Definition. A vertex x S T is critical for a chip configuration u if x ^ s and 

(7) u(x) < S C(x) | y is critical}. 

Proposition 3.1. A stable configuration u 6 SP(T) is recurrent if and only if 
equality holds in |?|) /or every critical vertex x. 
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A A A AA AA A 

Figure 2. A recurrent configuration on the wired ternary tree of 
height five. The top and bottom edges lead to the sink. Critical 
vertices are circled; if any of the circled vertices had fewer chips, 
the configuration would not be recurrent. 

Proof. If x is critical, then 

(8) u(x) + #{y S C(x) | y is not critical} < deg(x) — 1. 

Thus after chips are added as prescribed in the burning algorithm, inducting upward 
in decreasing distance to the root, if x ^ r is critical, its parent must topple before 
it does. In particular, if strict inequality holds in ([7]), and hence in pj, for some 
vertex x, that vertex will never topple, so u is not recurrent. 

Conversely, suppose equality holds in (J7]), hence in ([8]), for every critical x. Begin 
toppling vertices in order of decreasing distance from the root. Note that a non- 
critical vertex x satisfies 

(9) u(x) + #{?/ S C(x) | y is not critical} > deg(x). 

Inducting upward, every non-critical vertex topples once. Hence by equality in 
once all vertices other than the root are stable, every critical vertex x has cither 
toppled (if its parent toppled) or is left with exactly deg(a;) — 1 chips (if its parent 
did not topple). In particular, the root now topples, as it was given an extra chip 
in the beginning. Now if ir is a critical vertex that has not yet toppled, its parent 
is also such a vertex. Inducting downward from the root, since all of these vertices 
are primed with deg(x) — 1 chips, they each topple once, and u is recurrent. □ 

The principal branches of T are the subtrees T 1 , . . . ,T k rooted at the children 
r±, . . . ,rk of the root r of T. The wired tree T l includes an edge from to the 
sink; thus r 4 has the same degree in T l as in T, as the edge from to r has been 
replaced by an edge from to the sink. 
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If Ui is a chip configuration on T l , and a is a nonnegative integer, we will use the 



a 



notation for the configuration on T which has a chips at the root 

y m,...,u k J 

and coincides with itj on T l . The following result is an immediate consequence of 
Proposition |3.1| 



Lemma 3.2. Let u = ( 

V Ui,...,U k 

(i) If u is recurrent, then each Ui is recurrent. 

(ii) If U\, . . . , Uk are recurrent and a = k, then u is recurrent. 

Write S x for a single chip at a vertex x, and denote by x = (e + S x )° the recurrent 
form of 5 X . Note that if u is recurrent, then by §5§ 

(10) (u + x)° = (u + 6 x )°. 

Write (r) for the cyclic subgroup of S I P(T) generated by f , and {(ft, . . . , ft)) for 

the cyclic subgroup of j=1 5P(T i ) generated by the element (ri,...,ffe). As 
mentioned in the introduction, the following theorem can be expressed as the pair of 
short exact sequences (B, (EL with R(f) = (f) and R{f\. ■■,T^) = ((h, h))- 



The group H(T) appearing in both sequences is the quotient (111. 



Theorem 3.3. Let T 1 , . . . ,T k be the principal branches of T . Then 



k 



(11) SP(T)/{r) ~ SP(fi) /{(h, • • ■ , h)) 

i=l 

where r, are i/ie roots o/ T, T l respectively. 
Proof. Define : 5P(f ) 0* =1 5P(f l ) by 

< 12 > (....*.,«)»<«> 



Lemma 



3.2 



i) ensures this map is well-defined. Note that if is 

\vt,...,v k 

recurrent and 

ui,...,u k J \\ Vt,...,v k 

then by ( |10| , either b < k and itj = for all i; or b — k and the root topples, 
in which case = (u, + f^) for all i. Thus </> descends to a map of quotients 

: SP(T)/(r) - 0t x SP(fi)/{(h, . . . , f fc )). ^ 

By adding two configurations without allowing the root to topple, the configu- 
rations on each branch add independently, hence by (10 1 and Lemma 3.2 ii) 

-jf 

for some nonnegative integer j. Thus is a group homomorphism. Moreover 



a \ ( b \Y = (( k 
ut,...,u k J V vt, ■ ■ ■ ,v k J J \\ (u 1 + v 1 )°,...,(u k + v k y 



is surjective by Lemma 3.2 ii). Finally, to show injectivity, suppose that for some 
c > we have 

Ui = (vi + cfi)° , i=l,...,k. 
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Figure 3. Multiples (fcf )° of the root f in the wired ternary tree 
of height four. Each column vector represents a chip configuration 
which is constant on levels of the tree. 



Then from (|10| we obtain 
b 



Vl,...,Vk 

for some nonnegative integer d 



+ c(k + l)f 



a 

Ui,...,Uk 



df 



□ 



4. Regular Trees 



In this section we show that for regular trees, Theorem 3.3 can be strengthened 
to express SP{T) as a direct sum. 

Let T n be the regular tree of degree d and height n, and T n the graph formed 
by collapsing its leaves to a single sink vertex s, and adding an edge from the root 
to the sink. The chip configurations which are constant on the levels of T n form a 
subgroup of SP(T n ). If each vertex at height k has a& chips, we can represent the 
configuration as a vector (oi, . . . , a n _i). If such a recurrent configuration is zero on 
a level, all vertices between that level and the root are critical, so by Proposition [XT] 
they must have d — 1 chips each. The recurrent configurations constant on levels 
are thus in bijection with integer vectors (<zi, . . . , a n _i) with < at < d — 1, subject 
to the constraint that if dj = then ai = . . . = di—i = d — 1. 

The following lemma uses the lexicographic order on vectors given by a < b 
if for some k we have a„_i = b n —i, . . . , a^+i = &fc+i and ak < b^. In the cyclic 
lexicographic order on recurrent vectors we have also (d — 1, . . . , d — 1) < (d — 
l,...,d- 1,0). 

Lemma 4.1. Ifu,v are recurrent configurations on T n that are constant on levels, 
write u v if v immediately follows u in the cyclic lexicographic order on the set 
of recurrent vectors. Then for every integer k > 0, we have 

(kf)° ~* ((fc + l)f)°. 

Figure [4] demonstrates the lemma for a ternary tree of height 4. 

Proof. By ([To]) we have 

((k + l)f)° = (kr + d r )°. 

Thus if (kf)° = (ai, . . . , a„_i) with a\ < d—1, then ((fe+l)r)° = (ai+1, ci2, . . . , <2n-i) 
as desired. Otherwise, if not all equal d — 1, let j > 1 be such that ai = . . . = 
a,-_i = d — 1 and < d — 1. Adding a chip at the root initiates the toppling 
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\ fln-1 / \ a n -l J \ d n -l J \ d n -l / \ On-1 ) 

as desired. If all aj = d — 1 the cascade will travel all the way down, ending in 
(d — 1, . . . , d — 1, 0) as desired. □ 

Proposition 4.2. Let T n be the wired regular tree of degree d and height n, and 
let R{T n j be the subgroup of SP(T n ) generated by f. Then i?(T„) consists of all 
recurrent configurations that are constant on levels, and its order is 

(d- l) n - 1 



(13) 



d-2 



Proof. Since the identity element e is constant on levels, and the property of being 
constant on levels is preserved by stabilization, for any k > the configuration 

(kr)° = (e+k6 r )° 

is constant on levels. Conversely, by Lemma |4.1| any recurrent configuration con- 
stant on levels can be expressed as a multiple of f. The number of such configura- 



tions is the number of integer vectors of the form (d — 1, . . . , d — 1, 0, dj, . 
with 1 < j < n and 1 < di < d — 1 for each i = j, . . . , n — 1, which is 

n-l 



l). 



J=0 



' d-2 



Index the nonsink vertices of T n by words of length < n — 2 in the alphabet 
[d— 1] = {1, . . . , d— 1}. For i = 1, . . . , n — 2 let <Xj be the automorphism of T n given 
by 

ai{wi . . . w k ) = w 1 ...(w i + l)...w k 
with the sum taken mod d— 1; if k < i then <Ji(w) = w. Given a map a : [n — 2] — > 

[d — 1} let a a be the composition YYi=i a i ■ 

If a is an automorphism of the form cr Q , write era for the chip configuration 
au(x) = u(a~ 1 x). Given recurrent chip configurations u and v on T n , if xi, . . . , x m 
are the vertices that topple in the stabilization of u + v, then 

ra 

(u + v)° = u + v + 2J A^. . 

3=1 

Since crA x — A ax we obtain 

m 

a(u + v)° = au + av + A CTXj . 
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The configuration on the right side is stable, recurrent, and = era + era (mod A), 
so it is equal to (era + av)° . Thus er is an automorphism of the sandpile group. 

Proposition 4.3. Let T n be the wired regular tree of degree d and height n, and 
let R(T n j — (f) be the subgroup of SP[T n j generated by the root. Then 

with d—1 summands of SP(T n —i) on the right side. 
Proof. Define p : SP(T n ) -> SP(f n ) by 

(14) p( u )=Ud-l) 2 Yl 

\ a:[n-2]-*[d-l\ 

By construction, p(u) is constant on levels, so the image of p lies in R[T n j by 



Proposition 4.2 Given u E R(T n j, since u is constant on levels we have a a i 
for all a. Since there are (d— 1)"~ 2 terms in the sum (14 1, we obtain 

p(u) = ((d-l) n u)° = u, 

where the second inequality follows from (|1 3|. Thus i?(T n ) is a summand of 
SP(T n ), and the result follows from Theorem 13.31 □ 



Proposition 4.3 can fail for non-regular trees. For example, if T is the tree 
consisting of a root with two children each of which have three children (Figure 4) , 

then f = ^ ^ has order 10 and the element x — ^ ^ satisfies Ax = f, so 

x has order 40. The total number of recurrent configurations is 4 • 4 • 3 — 8 = 40, so 
SP(f) ~ Z/40Z, and R(f) ~ Z/10Z is not a summand. 

We now turn to the proof of Theorem |1.2| which can be written as 

sp(f n ) ~ z qn e 7L\z\ © z^; 1 ) © ... © z^"" 3 ^-^, 

where q n = 1 + a + . . . + a" -1 . 
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Proof of Theorem \l.S\ By Proposition |4.3| we have 

(15) SP(f n ) ~H n ® R(f n ) 
where 

(16) H n = SP{f n _ 1 )® a /D n _ 1 
and D n _i is the diagonal copy of i?(T„_i). 



By Proposition 4.2 we have R[T n ) — %*q n , so it remains to show that 

(17) H n ~ v i q -\ ffi z a q [ a -v e . . . e z^ 3 ^- 1 ). 

We will show this by induction on n; the base case n = 2 is trivial. Substituting 
( 15 1 into ( 16 1 gives 



tp 



By Proposition 4.2 we have i?(T„_i) ~ Z<?„-ii an d (17) now follows from the 
inductive hypothesis. □ 

5. Proof of Toumpakari's Conjecture 

As before, write a = d — 1 and 

q n = l + a+ .. . + a n ~ 1 . 

If p is a prime not dividing d(d — 1), let £ p be the least positive n for which p\q n . 
Then 

I p, if a = 1 (mod p) 

lordp(a), else. 

Here ord p (a) is the least positive k for which p|a fc — 1. Note that p\q n if and only 
if t p \n. The following result was conjectured by E. Toumpakari in 12J (where the 
factor of d — 2 was left out, presumably an oversight). 

Theorem 5.1. Let B n be the ball of radius n+ 1 in the d-regular tree, with leaves 
collapsed to a single sink vertex, but with no edge connecting the root to the sink. 
Let p be a prime not dividing d(d~ I), and let S p (n) be the Sylow-p subgroup of the 
sandpile group SP(B n ). Then 

(d(d-2)J2 m<n {d-l) m , if -1 (mod t p ); 

rankCS (n)) = < m ="( mo(li p) 

y py " \d{d-2)J2 m<n (d-l) m + d-l, ifn = -l(modtp). 

V m=n (mod t p ) 

Proof. Since each of the d principal branches of B n is isomorphic to T n+ i, by 
Theorem 13.31 we have 

SP(B n )/{f) ~ SP(f n _ +1 )®...(BSP(f n+1 ) 
(R(T, l+ i), . . . , i?(T„ + i)) 



4.2 



we have R(T n+ ij ~ Z 9n+1 , so from Theo- 



with d summands. By Proposition 
rem 11.21 

(18) SP(B n )/(f) ~ Z a qn+1 © z q a - 1){a+1) © z[ a -_\ )a(a+1) © ■ ■ • © Z^" 1)a " _2(a+1) . 

By Proposition 7.2 of [T2], the root subgroup (f) of 5P(B n ) has order d(d — l) n . 
Thus for p not dividing d(d ~ 1) the Sylow p-subgroup of SP(B n ) has the same 
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rank as that of the quotient SP(B n ) /(f). Each summand 1 qk in ( 18 ) contributes 1 
to the rank of S p (n) if t p \k and otherwise. If n ^ —1 (mod t p ), the total rank is 
therefore 

rank(5 p (n))= ^ {a - l)a n - k (a+ 1) 

2<k<n 

= d(d-2) J] (rf- 1 )" 1 - 

0<m<n-2 
m=n (mod t p ) 



In the case that n= — 1 (mod t p ), the first summand in (18) contributes an 

additional rank a = d — 1 to S p (n). □ 
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